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Abstract 

^vj , We consider a two component hybrid inflation model, in which two 

^ ' fields drive inflation. Our results show that this model generates an ob- 

servable non-gaussian contribution to the curvature spectrum, within the 
limits allowed by the recent WMAP year 3 data. We show that if one 
^*G ' field has a mass ri^ < 0, and an initial value (f> < O.OGMpi while the other 

\l I field has a mass rjo- > 0, and initial field value 0.5Mpi < a < Mpi then the 

non-gaussianity is observable with 1 < /nl < 1.5, but that /nl becomes 
much less than the observable limit should we take both masses to have 
the same sign, or if we loosened the constraints on the initial field values. 



o 
o 



c^ 



L ; 1 Introduction 

H , 

"^ I We consider hybrid inflation driven by two scalar fields and a, with the po- 

C^ . tential: 

-^; y = Fo(l + f0^ + f.^) (1) 

;_j ' where as usual the vacuum expectation value of the potential Vq dominates, 

770, r^cr are the masses-squared of the (f> and a fields respectively and may have 
either sign. We take Mpi — 1 throughout. This is similar to the case considered 
in Ref.|PI2], with the cr = trajectory, and since the corresponding |/nl| was 
found to be much less than 1, we aim to derive the non-gaussianity generated 
by a more general trajectory. We calculate the curvature perturbation and the 
resulting non-gaussianity using the AA^ formalism [2j (which is equivalent to 
using second-order cosmological perturbation theory as described in 131), for 
which: 

C(i,x) = 5N 

= NM^ + lN^JS<j),S(bJ (2) 
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Here N is the unperturbed number of e— folds from an initial time with field 
values (jji to a final time with energy density p, Ni = |^ and Nij = q^.^i. ■ 

We take the final time to be just before the end of inflation ^, and assume slow 
roll so that p = V. 

We measure the non-gaussianity via the amplitude of the bi-spectrum: 

where hi(kaH) ^ 1, k^^ is the scale under consideration, a is the scale 
factor, H = d/a is the Hubble parameter, and V(^ is the spectrum of curvature 
perturbations. This expression assumes that the S(j)i are gaussian and applies 
even if |/nl| ^ 1 "^171, the condition for /nl to eventually be observed. 



2 The Derivatives 

The amplitude of curvature perturbations defined in Eq. Q is related to the 
amount of expansion that occurs between an initial time, corresponding to a 
flat slice of space-time and a final time, corresponding to a space time slice of 
uniform energy density. Following the procedure of PP (which is equivalent to 
that used in Rcf. 7 ), the respective variations of the fields (p and a are defined 
purely by their initial values (denoted here as simply (/) and a) and the amount 
of inflation. 

We use the slow roll equation: 

3H(pn H TTi = (4) 

d(t)n 

where 0„ = d(j}n/dt. Rearranging this equation, and integrating over the 
period of inflation we get: 



" |§<i» (6) 

where we have used dN/dt ~ H (valid for slow roll) , and the * denotes the time 
when cosmological scales left the horizon. 

Recalling that 3H^ = ^n^i't'n) — Vq we have 4'{N) = 0exp(— r^^A^) and 
a{N) = aexp{—r]aN). Substituting these equations into Eq. |^ we have: 



V 



Vo (l - f 0' exp(-2r;^7V) - ^a^ exp(-2r;.7V)) (7) 



^by specifying a mechanism to end inflation, we could also calculate the contribution from 
ending as described in |5], but we do not consider this case here. 



differentiating Eq. iQ with respect to (j), while recalhng that Vi — and 
rearranging: 



dN 



m, 



?70(/)exp(-2iV?70) 



where 



p{(j), a, N) = r^lcf)'' exp(-27V77^) + r/^^ exp{-2Nri^) 
By differentiating Eq. ^ with respect to a we get: 

77<jCrexp(-2A^77c,) 



Defining: 



we get: 



N„ = 



(3 



3_2 . 



7(0, a, N) = 77^0" exp(-2iV77^) + 77>" exp(-2Ar77,) 



(8) 

(9) 
(10) 

(11) 



770 exp(-2iV77^) 773 02 exp(-47V7?0) 27 2 ,2 / . .r ^ ^1o^ 
iV00 = -^ '- ^ - 4^ -^ + -j^vlr exp{-ANr,^) (12) 



77^exp(-2JVr/^) ri^a^ expj-ANij^) 2j ^ 2 / aat \ nQ^ 

Naa = ^ 4 -^ + —T]^a expi-AN-q^) (13) 

and 

iV,. = N^, = 277,077.. exp(-^2iV(7y, + ,.)) ^ _ (^^ ^ ,^^)^) (^^^ 



3 Curvature Perturbation 

Substituting the equations from Scction[21into Eq. ^ we find that the curvature 
perturbation is given by the equation below. At first order, ^ is separable in 
terms of the individual field perturbations 6(f> and Sa, and involves a cross term 
{5(p6(j) at second order as well as the pure ((50^)^. 



c = 



2P 



27700exp(— 2A^77,)(50 + 277.0- exp(— 27V 7/cr)5cr 



3 J,3 



770 exp(-2Ar,y^) - 4^ exp(-4iV770) + S^t?^'/'' exp(-47V770) 



7730-3 



"f „2 2, 



77. exp(27V77.) - 4-^ exp(-47V77.) + 2-^77^.^ exp(-4iV77.) 



/3 



P' 



{5af 



2 ^ (^^ - iv, + V.) 



64)S<T 



(15) 



/] 
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Substituting Eqs.JSJ to (|14|l into Eq. (j2l we find that the non-gaussian contri- 
bution to the curvature spectrum is given by: 



/NL 



2^2 



exp(— 4A^?70) + ry^CT^ exp(— 4iVry(j) 



xl -(3 [~r]la^ exp(-6iV?7^) - Tyj^^^ exp{-6N7]^)] 

+2 [~2ria^ exp{-8mj„) - 2^^^ exp{~8m^) 
-2(77^ + V.)vlvl<l)^<y^ exp(-47V(?7, + r;^))] 
+2^ [r74a4exp(-8iV77,) +774^4 gxp(-87V770) 

+2772^2^^2 exp(-47V(77, + 77^))] I + • ■ • 



(16) 



We found tliat the second term is negUgible with respect to the first, so we 
do not inchidc it here. 



5 The Positive-Negative Mass Combination 

In this section we focus on taking 770- > while 77^ < 0. In this case the a field 
is pushing the inflaton towards the origin while the cf) field is pulling it away. 

To first order in field perturbations, wc find that the curvature perturbation 
Eq. (0) is dominated by the fluctuations in the negative mass field, since the 
fluctuations Sa are exponentially damped. 



■n<t>(p 



-2^ exp(2A^(7/0 - ■q„))5a 



(17) 



This argument could be used to simplify the non-gaussian term (|16|l with 
only a small loss of precision, especially for the case where 77^ ~ 1, and would 
reproduce the /nl for a single held model. However this is not the case for 
77o- -C 1, as then the exponent 2Nria ~ 1, and it becomes interesting to consider 
the affect of the Sa fluctuations on observation. 



"-i=2 j,. ..J ,,j: , :,_,_) ,,,^; (18) 



6 The Spectral Index 

For a multi-field inflaton the dominant term in the spectral index is defined by 



so by substituting Eqs. © and (^UJ), and using Eq. (QJ to calculate Vab we 
get: 

-3^.2 exp(-4iVr/0) + rjla^ exp(-47V77„) 
'77^02 exp(~4A^r/0) + ■qla'^ exp{~ANria) 

From the recent WMAP year three data 9 the range of allowed n — 1 for a 
negligible tensor fraction r at the 2a limit is: 

- 0.083 < 71 - 1< -0.02 (19) 

If we want the spectral index to fall within observational limits regardless 
of initial conditions, then we require the spectral index given in Eq. (|18(l to 
fall within the observational bounds Eq. H19|l independently of the initial field 
values. For the case where both rja^rj^ < then since Eq. (|18|l is a weighted 
sum of 2ri.i we require: 

h\<OMl (20) 

However, for the case where ry^ < and t]^ > 0, then as long as rj^ satisfies 
the limits set by Eq. H20|) then rja- can take on any value less than 1, since 
for large values of t]^ terms within the equations including it are exponentially 
damped. Also note that for positive field masses, n — 1 > which is ruled out 
by observation. 

7 Results for the Non-Gaussianity 

We begin by analyzing Eq. (|16|l for A^ = 50. Keeping within the range of 
rji allowed by Eq. 119|l . we have analysed the potential for initial field values 
ranging between (0 — 1), and values of the masses ranging between (—0.04 — 0) 
for 7^0 and (—0.04 — 0.15) for rj^. We then extracted the maximum value of /nl 
for each combination of masses and plot the results in Fi^ We repeat this 
procedure and extract ||/NL|ma£c for each combination of initial field values and 
plot these results in Fi^ We have also provided illustrations of the potential 
for the two scenarios considered in this paper in Fig^and^ For the case where 
77o- ^ and —0.04 < 77^ < 0, Eq. p6|l reproduces the results for a single field 
scenario. 

It is a far from easy task to extract any meaningful information from a four 
parameter model such as this one. However, both Figs. n]and[21are results of a 



full parameter space calculation of /nl, but FigQ] shows the mass dependence 
of |/nl|, while Fig|21 shows the initial field value dependance of |/nl| and by 
analyzing the figures we can see that for large a, small 0, and r^^ > 0,77^ < 
we get 1 < /nl < 1-5. Yet, if we were to consider small values of the fields 
(/), cr < 0.1, we see that /nl ^ 1 regardless of the values of the masses. 
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Figure 1: Contour plot of the absolute value of the maximum non-gaussianity 
generated for various combinations of masses, with initial conditions rang- 
ing between zero and the Planck mass. |/nl| ^ 1 corresponds to the four 
inner contours and the region outside the outermost contour corresponds to 
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<0.2. 
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Figure 2: Contour plot of the absolute value of the maximuni non-gaussianity 
generated for various combinations of initial field values, with masses ry^ ranging 
between (—0.04 — 0) and r/c ranging between (—0.04 — 0.15). The embedded 
plot is a magnification of the region in which |I/nl| > 0.2 appears, note that 
this corresponds to </> ^ 1 and 0.2 < a < 1. 




Figure 3: The potential Eq. (^ for rj^^rjcr > 0, in which the inflaton starts at 
a maximum and is pushed away from the origin by both fields. The solid black 
line corresponds to the unperturbed case a = 




Figure 4: The potential Eq. (^ for 77^ > and ri^ < 0, in which the (j) field 
pushes the inflaton away from the origin while the a field pulls it towards it. 
This case produces the more significant (yet still undetectable) values of /nl 
which appear in Figs^yil The solid black line corresponds to the unperturbed 
case cr = 

8 Discussion 

We summarise the main results in the table below, and by Consistent with n—1? 
we mean does it satisfy the limits on the spectral index in Eq. (|19|l . 



V<T 


Vcp 


a 





I/nlI 


Consistent 

with 

n-1? 


< r/CT < 1 


< 


< 1 


<1 


1 < I/nlI < 1.5 


yes 


< 1 


< 


any value 


any value 


single field case 
/nlI « 1 


yes 


<0 


<0 


any value 


any value 


/nlI « 1 


yes 


any value 


any value 


<1 


<1 


/nl| < 1 


yes 


> 


>0 


any value 


any value 


not 
calculated 


no 
n-1 > 


^0 


<0 


any value 


any value 


single field case 
I/nlI < 1 


yes 



For the first result listed r]cr -^ I, considering that an unperturbed trajectory 
corresponds to ct = (solid black lines in Figs and 0}, the non-gaussian 
contribution to the curvature spectrum is due to the slight deviation in the 



inflaton trajectory caused by the a field, and since —2Nr]a ^ 1, tiie effect 
'survives' inflation, leaving it's imprint on the spectrum. In the rj^ < I case, the 
roll to the minimum occurs rapidly at the beginning of inflation (by rapidly we 
mean on a time scale much smaller than the time scale of inflation), and thus 
the effect is negligible, and the results for the non-gaussianity are equivalent to 
a single field case |/nl — V<t>/'^- 

It is only fair to point out that we have not provided a physical motivation 
for this type of model with field values of order one, and since much effort was 
put into justifying multi- field chaotic models (as summarized in |EI|)i it is not 
clear whether we can write down our potential without justification. 
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Abstract 

^vj , We consider a two component hybrid inflation model, in which two 

^ ' fields drive inflation. Our results show that this model generates an ob- 

servable non-gaussian contribution to the curvature spectrum, within the 
limits allowed by the recent WMAP year 3 data. We show that if one 
^*G ' field has a mass ri^ < 0, and an initial value (f> < O.OGMpi while the other 

\l I field has a mass rjo- > 0, and initial field value 0.5Mpi < a < Mpi then the 

non-gaussianity is observable with 1 < /nl < 1.5, but that /nl becomes 
much less than the observable limit should we take both masses to have 
the same sign, or if we loosened the constraints on the initial field values. 



o 
o 



c^ 



L ; 1 Introduction 

H , 

"^ I We consider hybrid inflation driven by two scalar fields and a, with the po- 

C^ . tential: 

-^; y = Fo(l + f0^ + f.^) (1) 

;_j ' where as usual the vacuum expectation value of the potential Vq dominates, 

770, r^cr are the masses-squared of the (f> and a fields respectively and may have 
either sign. We take Mpi — 1 throughout. This is similar to the case considered 
in Ref.|PI2], with the cr = trajectory, and since the corresponding |/nl| was 
found to be much less than 1, we aim to derive the non-gaussianity generated 
by a more general trajectory. We calculate the curvature perturbation and the 
resulting non-gaussianity using the AA^ formalism [2j (which is equivalent to 
using second-order cosmological perturbation theory as described in 131), for 
which: 

C(i,x) = 5N 

= NM^ + lN^JS<j),S(bJ (2) 
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Here N is the unperturbed number of e— folds from an initial time with field 
values (jji to a final time with energy density p, Ni = |^ and Nij = q^.^i. ■ 

We take the final time to be just before the end of inflation ^, and assume slow 
roll so that p = V. 

We measure the non-gaussianity via the amplitude of the bi-spectrum: 

where hi(kaH) ^ 1, k^^ is the scale under consideration, a is the scale 
factor, H = d/a is the Hubble parameter, and V(^ is the spectrum of curvature 
perturbations. This expression assumes that the S(j)i are gaussian and applies 
even if |/nl| ^ 1 "^171, the condition for /nl to eventually be observed. 



2 The Derivatives 

The amplitude of curvature perturbations defined in Eq. Q is related to the 
amount of expansion that occurs between an initial time, corresponding to a 
flat slice of space-time and a final time, corresponding to a space time slice of 
uniform energy density. Following the procedure of PP (which is equivalent to 
that used in Rcf. 7 ), the respective variations of the fields (p and a are defined 
purely by their initial values (denoted here as simply (/) and a) and the amount 
of inflation. 

We use the slow roll equation: 

3H(pn H TTi = (4) 

d(t)n 

where 0„ = d(j}n/dt. Rearranging this equation, and integrating over the 
period of inflation we get: 



" |§<i» (6) 

where we have used dN/dt ~ H (valid for slow roll) , and the * denotes the time 
when cosmological scales left the horizon. 

Recalling that 3H^ = ^n^i't'n) — Vq we have 4'{N) = 0exp(— r^^A^) and 
a{N) = aexp{—r]aN). Substituting these equations into Eq. |^ we have: 



V 



Vo (l - f 0' exp(-2r;^7V) - ^a^ exp(-2r;.7V)) (7) 



^by specifying a mechanism to end inflation, we could also calculate the contribution from 
ending as described in |5], but we do not consider this case here. 



differentiating Eq. iQ with respect to (j), while recalhng that Vi — and 
rearranging: 



dN 



m, 



?70(/)exp(-2iV?70) 



where 



p{(j), a, N) = r^lcf)'' exp(-27V77^) + r/^^ exp{-2Nri^) 
By differentiating Eq. ^ with respect to a we get: 

77<jCrexp(-2A^77c,) 



Defining: 



we get: 



N„ = 



(3 



3_2 . 



7(0, a, N) = 77^0" exp(-2iV77^) + 77>" exp(-2Ar77,) 



(8) 

(9) 
(10) 

(11) 



770 exp(-2iV77^) 773 02 exp(-47V7?0) 27 2 ,2 / . .r ^ ^1o^ 
iV00 = -^ '- ^ - 4^ -^ + -j^vlr exp{-ANr,^) (12) 



77^exp(-2JVr/^) ri^a^ expj-ANij^) 2j ^ 2 / aat \ nQ^ 

Naa = ^ 4 -^ + —T]^a expi-AN-q^) (13) 

and 

iV,. = N^, = 277,077.. exp(-^2iV(7y, + ,.)) ^ _ (^^ ^ ,^^)^) (^^^ 



3 Curvature Perturbation 

Substituting the equations from Scction[21into Eq. ^ we find that the curvature 
perturbation is given by the equation below. At first order, ^ is separable in 
terms of the individual field perturbations 6(f> and Sa, and involves a cross term 
{5(p6(j) at second order as well as the pure ((50^)^. 



c = 



2P 



27700exp(— 2A^77,)(50 + 277.0- exp(— 27V 7/cr)5cr 



3 J,3 



770 exp(-2Ar,y^) - 4^ exp(-4iV770) + S^t?^'/'' exp(-47V770) 



7730-3 



"f „2 2, 



77. exp(27V77.) - 4-^ exp(-47V77.) + 2-^77^.^ exp(-4iV77.) 



/3 



P' 



{5af 



2 ^ (^^ - iv, + V.) 



64)S<T 



(15) 



/] 
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Substituting Eqs.JSJ to (|14|l into Eq. (j2l we find that the non-gaussian contri- 
bution to the curvature spectrum is given by: 



/NL 



2^2 



exp(— 4A^?70) + ry^CT^ exp(— 4iVry(j) 



xl -(3 [~r]la^ exp(-6iV?7^) - Tyj^^^ exp{-6N7]^)] 

+2 [~2ria^ exp{-8mj„) - 2^^^ exp{~8m^) 
-2(77^ + V.)vlvl<l)^<y^ exp(-47V(?7, + r;^))] 
+2^ [r74a4exp(-8iV77,) +774^4 gxp(-87V770) 

+2772^2^^2 exp(-47V(77, + 77^))] I + • ■ • 



(16) 



We found tliat the second term is negUgible with respect to the first, so we 
do not inchidc it here. 



5 The Positive-Negative Mass Combination 

In this section we focus on taking 770- < while 77^ < 0. In this case the a field 
is pushing the inflaton towards the origin while the cf) field is pulling it away. 

To first order in field perturbations, wc find that the curvature perturbation 
Eq. (0) is dominated by the fluctuations in the negative mass field, since the 
fluctuations Sa are exponentially damped. 



■n<t>(p 



-2^ exp(2A^(7/0 - ■q„))5a 



(17) 



This argument could be used to simplify the non-gaussian term (|16|l with 
only a small loss of precision, especially for the case where 77^ ~ 1, and would 
reproduce the /nl for a single held model. However this is not the case for 
77o- -C 1, as then the exponent 2Nria ~ 1, and it becomes interesting to consider 
the affect of the Sa fluctuations on observation. 



"-i=2 j,. ..J ,,j: , :,_,_) ,,,^; (18) 



6 The Spectral Index 

For a multi-field inflaton the dominant term in the spectral index is defined by 



so by substituting Eqs. © and (^UJ), and using Eq. (QJ to calculate Vab we 
get: 

-3^.2 exp(-4iVr/0) + rjla^ exp(-47V77„) 
'77^02 exp(~4A^r/0) + ■qla'^ exp{~ANria) 

From the recent WMAP year three data 9 the range of allowed n — 1 for a 
negligible tensor fraction r at the 2a limit is: 

- 0.083 < 71 - 1< -0.02 (19) 

If we want the spectral index to fall within observational limits regardless 
of initial conditions, then we require the spectral index given in Eq. (|18(l to 
fall within the observational bounds Eq. H19|l independently of the initial field 
values. For the case where both rja^rj^ < then since Eq. (|18|l is a weighted 
sum of 2ri.i we require: 

h\<OMl (20) 

However, for the case where ry^ < and t]^ > 0, then as long as rj^ satisfies 
the limits set by Eq. H20|) then rja- can take on any value less than 1, since 
for large values of t]^ terms within the equations including it are exponentially 
damped. Also note that for positive field masses, n — 1 > which is ruled out 
by observation. 

7 Results for the Non-Gaussianity 

We begin by analyzing Eq. (|16|l for A^ = 50. Keeping within the range of 
rji allowed by Eq. 119|l . we have analysed the potential for initial field values 
ranging between (0 — 1), and values of the masses ranging between (—0.04 — 0) 
for 7^0 and (—0.04 — 0.15) for rj^. We then extracted the maximum value of /nl 
for each combination of masses and plot the results in Fi^ We repeat this 
procedure and extract ||/NL|ma£c for each combination of initial field values and 
plot these results in Fi^ We have also provided illustrations of the potential 
for the two scenarios considered in this paper in Fig^and^ For the case where 
77o- ^ and —0.04 < 77^ < 0, Eq. p6|l reproduces the results for a single field 
scenario. 

It is a far from easy task to extract any meaningful information from a four 
parameter model such as this one. However, both Figs. n]and[21are results of a 



full parameter space calculation of /nl, but FigQ] shows the mass dependence 
of |/nl|, while Fig|21 shows the initial field value dependance of |/nl| and by 
analyzing the figures we can see that for large a, small 0, and r^^ > 0,77^ < 
we get 1 < /nl < 1-5. Yet, if we were to consider small values of the fields 
(/), cr < 0.1, we see that /nl ^ 1 regardless of the values of the masses. 
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Figure 1: Contour plot of the absolute value of the maximum non-gaussianity 
generated for various combinations of masses, with initial conditions rang- 
ing between zero and the Planck mass. |/nl| ^ 1 corresponds to the four 
inner contours and the region outside the outermost contour corresponds to 
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Figure 2: Contour plot of the absolute value of the maximuni non-gaussianity 
generated for various combinations of initial field values, with masses ry^ ranging 
between (—0.04 — 0) and r/c ranging between (—0.04 — 0.15). The embedded 
plot is a magnification of the region in which |I/nl| > 0.2 appears, note that 
this corresponds to </> ^ 1 and 0.2 < a < 1. 




Figure 3: The potential Eq. (^ for rj^^rjcr > 0, in which the inflaton starts at 
a maximum and is pushed away from the origin by both fields. The solid black 
line corresponds to the unperturbed case a = 




Figure 4: The potential Eq. (^ for 77^ > and ri^ < 0, in which the (j) field 
pushes the inflaton away from the origin while the a field pulls it towards it. 
This case produces the more significant (yet still undetectable) values of /nl 
which appear in Figs^yil The solid black line corresponds to the unperturbed 
case cr = 

8 Discussion 

We summarise the main results in the table below, and by Consistent with n—1? 
we mean does it satisfy the limits on the spectral index in Eq. (|19|l . 



V<T 


Vcp 


a 





I/nlI 


Consistent 

with 

n-1? 


< r/CT < 1 


< 


< 1 


<1 


1 < I/nlI < 1.5 


yes 


< 1 


< 


any value 


any value 


single field case 
/nlI « 1 


yes 


<0 


<0 


any value 


any value 


/nlI « 1 


yes 


any value 


any value 


<1 


<1 


/nl| < 1 


yes 


> 


>0 


any value 


any value 


not 
calculated 


no 
n-1 > 


^0 


<0 


any value 


any value 


single field case 
I/nlI < 1 


yes 



For the first result listed r]cr -^ I, considering that an unperturbed trajectory 
corresponds to ct = (solid black lines in Figs and 0}, the non-gaussian 
contribution to the curvature spectrum is due to the slight deviation in the 



inflaton trajectory caused by the a field, and since —2Nr]a ^ 1, tiie effect 
'survives' inflation, leaving it's imprint on the spectrum. In the rj^ < I case, the 
roll to the minimum occurs rapidly at the beginning of inflation (by rapidly we 
mean on a time scale much smaller than the time scale of inflation), and thus 
the effect is negligible, and the results for the non-gaussianity are equivalent to 
a single field case |/nl — V<t>/'^- 

It is only fair to point out that we have not provided a physical motivation 
for this type of model with field values of order one, and since much effort was 
put into justifying multi- field chaotic models (as summarized in |EI|)i it is not 
clear whether we can write down our potential without justification. 
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